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, We present a numerical study of sidebranching of a solidifying dendrite by means of a phase- 

field model. Special attention is paid to the regions far from the tip of the dendrite, where linear 
theories are no longer valid. Two regions have been distinguished outside the linear region: a 
first one in which sidebranching is in a competition process and a second one further down where 
branches behave as independent of each other. The shape of the dendrite and integral parameters 
CO ' characterizing the whole dendrite (contour length and area of the dendrite) have been computed 

, and related to the characteristic tip radius for both surface tension and kinetic dominated dendrites. 

Conclusions about the different behaviors observed and comparison with available experiments and 
theoretical predictions are presented. 

o . 

, PACS numbers: 81.10.-h, 81.10.Aj, 81.10.Dn, 64.70.Dv 

I. INTRODUCTION 

' The generation of dendritic patterns arises in different noncquilibrium situations. [1, 2, 3, 4, 5, 6, 7, 8]. The case of 
dendrites appearing during the solidification of a melt has long provided an archetypical example of pattern forming 
I . system, in which the underlying physics is well known. Nonetheless, from the theoretical point of view it has posed 
a number of nontrivial questions on the selection of the final growth mode, which in a large amount has driven the 
research on the effects of nonlinearities, anisotropics and fluctuations in interfacial pattern formation. These questions 
^ also have an applied interest, since solidification is one of the most common methods to produce materials. It is well 
known that the details of the dendritic pattern (and in particular its associated scales) appearing during growth, 
determine the microstructure of the grown solid, which in turn is responsible to a large degree of its final (mechanical 
and electrical) properties [9]. 

. In this context an increasing attention is focussed on the shape of the growing dendrite and on sidebranching 
^ ; [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37], which 
\^ . corresponds to the apparition and growth of secondary branches at both sides of the dendrite. Sidebranching activity 
0^ shows different behaviors depending on the distance to the tip of the dendrite [16, 34, 35]. In the zone closer to the 
tip sidebranches are born as a convective instability of the dendrite and grow linearly. Further down from the tip, 
sidebranches are usually much more developed and a competition process between branches takes place mediated by 
the interaction between their diffused fields. Much further from the tip, the competition has finished and winner 
I branches grow as free dendrites while the growth of looser branches is inhibited, 
g ; Many theoretical [12, 13, 14, 15, 16, 17, 18, 19, 20] and experimental [21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 
i' 33, 34, 35] studies of the region close to the tip have been carried out in the recent years. A common point in the 
'■^ study of the linear region has been the characterization of sidebranching by means of its amplitude and wavelength. 
^ , The growth of sidebranches in the regions further down from the tip presents a behavior very different from that 
O • of the linear regime. One finds first a region where branches compete interacting through the expelled heat. This 
[ gives rise to an irregular growth of sidebranches which is difficult to characterize by their amplitude and wavelength, 
^ 1 since these quantities are not longer well defined outside the linear region. Experiments carried out with different 
■ substances [25, 26, 30] have shown that sidebranching in this region is self-similar and that geometrical parameters 
, can be scaled by the tip radius R. This nonlinear region and its associated self-similar growth have their limits at 
d ' distances to the tip of the order of the diffusion length, i.e. z/R « 1/Pe [16], where Pe is the Peclet number. 
Further down, sidebranches behave like dendrites themselves. 

The nonlinear region has been intensively studied in experiments with xenon dendrites by Hiirlimann et al. [26] and 
with succinonitrile dendrites by Li and Beckermann [30, 31]. In particular, the shape of the sidebranching envelope 
was studied in Ref. [30] by measuring the distance X from the axis of the dendrite to the tip of active sidebranches 
(defined as those branches longer than all the other branches closer to the tip) versus the distance Z to the tip along 
the axis of the dendrite (see Fig. 1). Values of X and Z were computed from the image of the dendrite projected on 
a plane, and far from the tip the relation X/R = 0.668(Z/i?)''-^^^ was obtained. 

It was proposed in Refs. [26, 30] an alternative set of integral parameters in order to describe the complex shape of 
a dendrite as a whole and the nonlinearities of dendritic solidification. Parameters characterizing independent parts 
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of the dendrite (e.g. amplitude and wavelength of the sidebranching) do not take into account the interaction of the 
sidebranches through the diffusion field. Nonlinear effects such as e.g. coarsening make unclear which sidebranches 
should be included in the measurement of the wavelength and which others should not. Instead, the contour length 
U, the projection area F and the volume of a dendrite appear to be more appropriate. 

It was found in the earlier experimental work of Refs. [26, 30] that the projection area varied linearly with the 
contour length and the corresponding slope M accomplished M/R = constant. Similar results were found in early 
simulations [36] , but this was shown to be an effect of reflecting boundary conditions strongly affecting the dendrite 
[37] . This suggests that the experimental observations could have been affected by the diffusion field of other close 
dendrites or other growing morphologies. In more recent experimental [31] and numerical [37] works it was obtained 
that F/ {UR) was not a constant in the nonlinear regime. 

The projection area showed two different behaviors (F/R^ = 0.847 (Z/Ry-^^^ for Z/R < 30 and F/R^ = 
0.578(Z/i?)^'^^ for Z/R > 30) in three-dimensional succinonitrile dendrites [31]. However, in two-dimensional am- 
monium bromide dendrites [22], the area F was found to vary over three orders of magnitude as Z^'^, as it would 
have happened if dendrites had had a smooth parabolic shape. As regards the variation of the contour length with 
the distance to the tip, only data corresponding to three-dimensional dendrites are available, where two behaviors are 
distinguished {U/R = 0.887(Z/i?)i "6 for Z/R < 20 and U/R = 0.378{Z/Ry-^° for Z/R > 40) [31]. 

An additional question is whether strong undercoolings can produce qualitative changes in sidebranching character- 
istics. It is well known that increasing undercooling the growth can switch from a regime dominated by surface tension 
to a regime dominated by kinetic effects. This was already predicted theoretically in Ref. [4]. When anisotropies 
of both effects favor different directions, changes in the growth directions of both dendrite and branches occur by 
changing undercooling. Even if these anisotropies are in the same directions, the behavior of the tip radius and 
velocity can present abrupt changes. A numerical evidence of such changes can be found in Ref. [45]. 

In this paper we present a study of sidebranching by means of a phase-field model for moving solid-liquid interfaces 
[38, 39, 40, 41, 42, 43, 44]. We consider sidebranching generated by selective amplification of fluctuations near the tip 
of a free growing dendrite. In particular, we focus on the nonlinear zone, including both the region where competition 
occurs and further down where sidebranches behave as free growing dendrites. Characterization is performed working 
out the shape of the dendrite by means of its envelope, and calculating the integral parameters. We have varied 
undercooling in a large range, in particular reaching relatively high values of the undercooling. This has permitted 
on the one hand, due to reduction in diffusion length, to access the region of free growing sidebranches far from the 
tip, and on the other hand to reach the kinetic regime of growth. 

This paper is organized as follows: In Sec. II we present the classical sharp interface model that characterizes a 
solidification system, the phase-field model and the numerical procedure used in this work. In Sec. Ill we present the 
results of simulations. We particularize in the effect of varying undercooling, and in the differences between different 
zones of the dendrite. Detailed characterization of the whole dendrite is performed by its shape and by computing 
the integral parameters. Finally, concluding remarks are presented in Sec. IV. 

II. MODEL AND NUMERICAL PROCEDURE 

The free solidification of a pure substance can be described by the sharp-interface model [1], which relies on 
the heat diffusion equation together with two boundary conditions at the interface, namely heat conservation and 
Gibbs-Thomson (local equilibrium) equation: 

dT 

— = DV'T, (2.1) 
ii;„ = £>Cp[(V„r)s-(V„T)i], (2.2) 



interface = Tm - -^H9) + a" {0)]n - VnPk{0). (2.3) 

In these equations T is the temperature (Tm being the melting one), D is the diffusion coefficient (D = k/cp, being 
k the heat conductivity and Cp the specific heat per unit volume), L is the latent heat per unit volume, u„ is the 
normal velocity of the interface, V„ is the normal derivative at the interface (S and L referring to solid and liquid 
respectively), (t{9) is the anisotropic surface tension (where 9 is the angle between the normal to the interface and some 
crystallographic axis) and k, is the local curvature of the interface. /3/c(6') is an anisotropic kinetic term, introduced 
into the Gibbs-Thomson Eq. 2.3 to account for a linear nonequilibrium correction. 
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The results of simulations presented below have been obtained by means of a phase-field model. These kind of 
models have received increased attention during the last years [44] . One of their main features is the introduction 
of an additional non-conserved scalar order parameter or phase-field (p, whose time evolution equation is coupled 
with the heat diffusion equation through a source term in order to take into account the boundary conditions at the 
interface. The phase-field takes constant values in each of the bulk phases (in our case, = in the solid and 0=1 
in the liquid) changing continuously between them over a transition layer, the interfacial thickness e. The equations of 
the model are then constructed in such a way that they converge to the sharp-interface dynamics of Eqs. (2.1,2.2,2.3) 
in the limit of vanishing e. Hence this parameter controls the convergence to the sharp-interface limit. 

The corresponding equations for the time evolution of the phase-field and the dimensionless temperature can be 
written in the following form [39]: 



— e 



dx 



dt 



dy 



1 



+ 30e/3AM^!>(l - cj)) 



d_ 
dy 



mv {0) 



dx 



+ e2V [r?'(e)V(/.] 



(2.4) 



^ + ^ (30</)2 - 60</.=^ + 30</)4) -^ = y\ + i;{x, y, t) (2.5) 

where u(r,t) is the diffusion field and A = CpAT/L is the dimensionless undercooling. Lengths are scaled by some 
arbitrary reference length u), while times are scaled by uP' jD. In these equations Q is the angle between the a;-axis 
and the gradient of the phase-field. t]{Q) = (T{9)/a{0) is the anisotropy of the surface tension. t{6) is given by 

^^77(0)/3/c(0), so the anisotropy of the kinetic term is given by T{9)/r](9). /3 is equal to and dg = CpTMcr{0)/L'^ 
is the capillary length. 

A source of fluctuations is introduced through the additive term tp in the heat equation. It was demonstrated [20] 
that sidebranching induced by this kind of noise qualitatively reproduce the characteristics of the (thermodynamical) 
internal noise, which makes it appropriate for the study of sidebranching. In our two-dimensional simulations the 
noise term is evaluated at each cell (i, j) of lateral size Ax as Ivij, where / denotes the amplitude of the noise, and 
rij is an uncorrelated uniform random number in the interval [—0.5,0.5]. The phase-field model equations have been 
solved on rectangular lattices using first-order finite diff^c^renccs on a uniform grid with mesh spacing Ax. An explicit 
time-differencing scheme has been used to solve the equation for cj), whereas for the u equation the alternating-direction 
implicit method was chosen [46]. The kinetic term has been taken as isotropic, which leads to t{6) = mr]{0) with 
constant m. A four-fold surface tension anisotropy r/(0) = 1 + 7cos(40) has been considered. 

The growth morphologies have been obtained by setting a small vertical seed (0 = 0, u = 0) in the center of the 
bottom side of the system and imposing 0=1 and w = — 1 on the rest of the system. Symmetric boundary conditions 
for and u have been used on the four sides of the system. Special care has been taken to employ large enough 
system sizes to avoid any influence of boundary conditions on the results presented along this paper. 

We have used a set of phase fleld model parameters that gives rise to a growing needle without sidebranching when 
no noise (/ = 0) is added to the simulations. This assures us that the sidebranching observed when / is not 
due to numerical noise. The fixed parameters for all the simulations have been f3 = 320, 7 = 0.045, m = 16 and 
e = 3.75 X 10"'^. The value of A has been varied in the range 0.44 — 0.65. The noise amplitude was kept constant 
(J = 16) in all the simulations and the time and spatial discretizations used were At = 1.25 x 10~^ and A.r = 0.0125. 

Under these conditions, the obtained morphologies have been dendrites with three main arms growing from the 
seed, one in the vertical {y) direction and two in the horizontal (x) one. We have focusscd on the sidebranches which 
grew perpendicular to the vertical arm. Thus, in order to get rid of the infiuence of the diffusion field of the horizontal 
arms on these sidebranches, we have been forced to run long simulations and only observe an area at a fixed distance 
to the tip. In order to avoid working with unnecessarily large systems, we have performed periodic shifts of the 
complete system, practice that has been checked it did not affect the results of the simulation. If Fig. 1 we show an 
example of a typical grown dendrite. 



III. RESULTS AND DISCUSSION 



A first series of simulations have been performed exploring the effect of undercooling on the tip radius. This has 
been measured by computing R = (j)x/4>yy at the tip of the dendrites [39]. The aim has been to identify different 
regimes of growth in a large range of A. In our case, anisotropy is only considered in surface tension. Thus, when 
the undercooling is changed no change in growth direction is expected although the behavior of the tip radius and 
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FIG. 1: Example of dendrite obtained at A = 0.575. Distances Z and X of active sidebranches used to characterize the shape 
of the dendrite, and the integral parameters contour length (?7) and surface area (F), are indicated. 
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FIG. 2: Tip radius of the dendrite vs. undercooling. 



velocity may vary. In particular in the surface tension controlled regime the tip radius should decrease (together with 

increase of tip velocity) by increasing undercooling. On the contrary in the kinetic regime, with isotropic kinetic term, 
one expects larger tip radius at higher velocities. In Fig. 2 it is shown the behavior of the tip radius as a function 
of the undercooling. It can be clearly observed a change of the behavior around the value 0.575. Thus, by choosing 
appropriate values of A we can select both regimes of growth. The existence of these two different regimes can be 
also confirmed by looking at the behavior of the tip velocity as a function of the undercooling or the Peclet number. 

We have looked at the shape of the studied dendrites by computing the coordinates (X, Z) , as defined above 
following Ref. [30]. Thus, only data of active branches are taken into account, i.e. branches longer than any other 
closer to the tip. 

Fig. 3 shows the plot of ZjR vs. X/R for three different values of the undercooling (A = 0.48,0.55,0.625). 
Representation for each A contains data from eight different times, which explains the slight dispersion of points. 
It can be distinguished in Fig. 3 the existence of two regimes for each undercooling, which is more evident as the 
undercooling is increased. For small values of Z/ R. similar behaviors are found for all the undercoolings. However, 
from a certain value of Z/R, X/R depends very much on A. The transition region between both regimes is not clear 
enough to permit the precise location of a crossover point in this figure. 

In order to better characterize the two observed regimes, it is shown in Fig. 4 the log- log plot of X/R vs. Z/R for 
the smallest and the largest undercoolings presented in Fig. 3. It can be observed from Fig. 4 a clear change in the 
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FIG. 3: Plot of ZjR vs. X/R for the sidebranches which are larger than any others closer to the tip. Symbols +, x and * 
correspond to A = 0.48, 0.55 and 0.625, respectively. 
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FIG. 4: Log-log plot of X/R vs. Z/R. for the sidebranches which are larger than any others closer to the tip. Symbols x and 
+ correspond to A = 0.48 and 0.65, respectively. In each case, regions A and B are indicated. 



behavior of X/R in the regions of Z/R around 80 and 40 for the case of small and large undercoolings, respectively. 
This suggests that these regions separate two zones A and B (see Fig. 4) where sidebranching is in different regimes. 
For A < 0.48 it is difficult to distinguish these two regimes because their slopes are very similar and it is not possible 
to determine a transition region. When the undercooling is increased it is found that the transition region is closer 
to the tip, which is consistent with the fact that at larger undercoolings a more developed sidebranching is obtained. 

The behavior of X/R in region A is not exactly the same for all the considered undercoolings and is given by a 
straight line in the log-log plot. By comparing data from different undercoolings, it can be observed that the set of 
points in region A lies at larger values of X/ R in the case of larger A. This is consistent with observations reported 
in Ref. [34, 35], where the exponents a calculated for each single branch in x ^ f^, x being branch length and t 
being time, were systematically smaller in branches grown in lower undercooling conditions. According to this, at any 
value of the undercooling and at small Z/ R, points in X/ R{Z/R) representation are less dispersed than at large Z/ R 
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because in the region closer to the tip it is still too soon to see the effects of the difference in the exponent a of each 
branch. Thus, the dispersion of points must increase with Z / R, as can also be observed in Fig. 3. 

The behavior of X/R in region B depends on A. When fitting the points of region B to X/R ~ {Z/R)"" we found 
that values of a tended to 1 for increasing A (in fact, from A > 0.525, variations in a are very small). In Fig. 5 it 
is shown a dendrite grown at A = 0.6 where regions A and B can be clearly distinguished. It can be observed that 
the angle 7 formed by the line joining the tips of active sidebranches and the axis of the main arm is smaller in the 
region A and, typically for this range of undercoolings, its value is very close to 45 degrees in the region B. In other 
words, in region B sidebranches grow at the same velocity that the main tip, i.e. grow as free dendrites. 
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FIG. 5: Dendrite grown at A = 0.6 where regions A and B are indicated. Angle 7 is always larger in regions B, where its 
value tends to 45 degrees as the undercooling is increased. 



The observation of these two regimes reveals a significant difference with experiments presented in Ref. [30] (sec 
Fig. 7 there), where only one regime was observed, and the angle formed by the axis of the main arm and the line 
joining the tips of the branches were always considerably smaller than 45 degrees. This has to be related to the small 
values of undercooling used in these experiments. The diffusion lengths associated with such a slow growth are very 
large, and even in the furthest from the tip region considered in the experiments the process of competition between 
branches was not finished yet. On the contrary, diffusion length in simulations is short due to the large undercoolings 
used. Branches can grow as free dendrites as long as distances between active sidebranches (which increase with the 
distance to the tip due to the competition process) are larger than the interaction scale given by the diffusion length 
(which is reduced for larger growth velocities). This results to be the condition for the zone B to be observed. 

We have also measured the integral parameters (contour length U and area F) of our two-dimensional dendrites. 
As shown in Fig. 1, f7 is the length of the contour of the dendrite measured from the tip to a distance Z along the 
axis, while F is one half of the area of the dendrite. Both magnitudes have been measiircd for different values of the 
dimensionless undercooling. Considering the origin of coordinates at the tip, the contour length and the area have 
been calculated from the coordinates of the dendrite contour by 

n 

U = ^[(Zi+1 - Zif + (Xi+1 - Xi)2]i/2 (3.1) 

i=l 

and 

F = ±^l2±l±lA^Z,^,-Z,), (3.2) 

where n corresponds to each distance to the tip for which we calculated U and F. The fact that the shape of 
sidebranches is rather irregular and that their growth is not always perpendicular to the y-axis, makes it difficult to 
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define U and F in a unique way everywhere as a function of Z . In order to better define both functions and following 
Ref. [30], we have only considered U and F for the values of Z corresponding to valleys between two neighboring 
sidebranches. 

It is shown in Fig. 6 the log-log plot of the normalized value of the surface area as a function of the normalized 
value of the distance to the tip of the dendrite at A = 0.525 with data taken at three different times. The same 
representation for the rest of the considered undercoolings shows the same behavior and only points far from the tip 
at larger A are slightly dispersed. 
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FIG. 6: Log-log plot of F/i?^ vs. Z/R, F being the surface area, for A = 0.525. 



F/B? vs. {Z/R) follows a power-law (a ~ 6°) where c is always around 1.5, although a slight tendency to increase 
with A is also observed. The value of c found in the simulations completely coincides with that found in Ref. [22] 
for the growth of ammonium bromide crystals in two dimensions. However, in Ref. [30, 31] the representation of 
F/B?{Z/R) showed two regimes of power-law behavior with different exponents. In principle, we should not expect 
to find the same exponent in our simulations taking into consideration that these experiments were three-dimensional 
and it was measured the projection area of the dendrite. In this case, the existence of different regimes near and far 
from the tip of the dendrite was attributed to the different effect of coarsening. This effect is also present in two 
dimensions, but the fact that only one regime is observed in the plot of F /B? makes us conclude that the manifestation 
of the coarsening effect is less dramatic in 2D that in 3D. As regards to the prefactor of the power-law fitting, it shows 
a similar behavior to that of the tip radius, that is, it decreases when the undercooling is increased up to A < 0.55 
and it increases for larger A. 

The behavior of the normalized contour length as a function of the distance to the tip for A = 0.44 and 0.6 is shown 
in Fig. 7. As it happened in the plot of the shape of the dendrite, the behavior changes after a transition region, 
being the variation of [//i? larger in the regions further down from the tip. 

In the region closer to the tip, it is found [//i? ~ {Z/Ry, which coincides with the behavior found in the linear 
regime of experiments in Ref. [31]. As one should expect, simulation results show that the linear region is larger for 
smaller undercoolings. This was not observed in Ref. [31], probably due to the employed range of undercoolings. 

After the linear region, there is a transition region which is followed by the nonlinear region, as it was observed in 
the experiments (see Fig. 5 in Ref. [31]). The transition region in the U / R{Z / R) plot is at smaller values of Z/R than 
the transition region in the X/R{Z/R) plot (Fig. 4) for all the considered undercoolings. In fact, the change in the 
behavior of the contour length and that of the shape of the dendrite (the envelop of it) provide different information of 
the sidebranching activity. In the case of U, the linear and nonlinear regimes are associated to low and high developed 
perturbations of the interface respectively. In the nonlinear region, both active and non-active branches contribute to 
the calculation of U. The enhanced growing of the active branches observed in the nonlinear region is accompanied by 
coarsening, process in which the shrinking of the shorter branches reduces the total increasing of the contour length. 
As a result any nonlinearity has effect on the behavior of U. On the contrary the shape X/R{Z/R) is calculated 
through the active branches only. X/R{Z/ R) is then associated to the effect on the winning branches of competition, 
and to these larger branches reaching colder regions and hence growing faster. As a result, dispersions of X/ R values 
are rather large, and the change of behavior more difficult to locate with a tendency to occur inside the nonlinear 
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FIG. 7: Log- log plot of U/R vs. Z/R, U being the contour length. Symbols x and + correspond to A = 0.44 and 0.6, 
respectively. Lines indicate the fitting of points in each region. 



region. The behavior of U/R in the nonlinear region can be fitted by a powder- law (see Fig. 7), although the values 
of the prefactor and the exponent depend on the undercooling. As regards to the prefactor, its variation with A is 
very similar to that of the tip radius, that is, it decreases up to A ^ 0.55 (surface tension dendrites) and it increases 
at larger undercoolings (kinetic dendrites). The largest value of the prefactor is 0.101 for A = 0.44, which is very far 
from the value obtained in Ref. [31]. The divergence is probably related to the different ranges of undercoolings used 
in simulations and experiments, but the difference in dimensions could also play a role. 

As regards to the exponent in the fitting of U/R, it increases with the undercooling from 1.57 to 1.89. The value 
for small A is very similar to the unique value (1.50) obtained in experiments [31]. Again, the fact of having found 
many exponents in the simulations and only one in the experiments could be associated to the different ranges of A 
used. It implies that difi^usion length considerably vary between simulations and experiments. The influence of the 
diffusion length on the competition process between branches that takes place in the nonlinear region determines the 
evolution of branches and consequently the behavior of the contour length. 

By combining results of the integral parameters in the linear regime, it is found F/ (U R) ~ {Z/ i?)"'^ , which coincides 
with the experiments [31]. The same exponent in the nonlinear regime varies from —0.07 to —0.39, differing very 
much from the experiments. Thus, the similarities between our results and the experimental ones in Ref. [31] remain 
mainly in the linear region and in the opposition to previous studies [26, 30], where F/UR was found to be constant 
in the nonlinear regime. 



We have presented a numerical study of the shape and sidebranching in regions at different distances from the tip 
of a solidifying dendrite by means of a phase-field model with a non-c:onserved noise term. We have characterized 
the dendrite by using the integral parameters and we have focussed in dendrites grown in both surface tension and 
kinetic regimes. 

The behavior of the shape of the dendrite has been foimd to depend on the undercooling in the considered range. 
The different diffusion lengths make the competition process between sidebranches to differ and thus the final shape 
of the dendrite is affected. The region where the competition process is taken place (A) and that where it is finished 
and sidebranches evolve like free dendrites (B) have been clearly distinguished in our simulation results. The behavior 
observed in B is in agreement with theoretical predictions [16]. On the other hand, the main divergence with the 
available experiments [30] is precisely the existence of two behaviors in the nonlinear region. This discrepancy may 
be explained by the different range of undercoolings considered in experiments and simulations. In our simulations 
undercooling (and hence tip velocity) is larger, so diffusion length is smaller and the transition between both zones 
is expected to occur closer to the tip, becoming observable. Note that additional increasings of undercooling would 
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reduce further the size of zone A, which then could not be considered as a separate scahng region. The area of the 
dendrite presented a unique behavior for aU the considered undercoohngs. As one would expect, it coincides with that 
of two-dimensional dendrites [22], although there is a slight discrepancy with three-dimensional dendrites, especially 
in regions far from the tip. The behavior F/B? ~ {Z/Ry -^ found in our simulations is the same as if the dendrites had 
a smooth parabolic shape. Thus, we can conclude that the area of two-dimensional dendrites is basically independent 
of the appearance and competition of sidebranches, and that situation does not depend of the diffusion length of the 
system. 

The behavior of the contour length presents two differentiated regimes in the linear and nonlinear regions. The 
exponents found in the linear region arc in agreement with experiments [31], while the discrepancies appeared in the 
nonlinear region could come from the range of undercoolings or the dimensionality. 

We have found that the behavior of the contour length changes in regions closer to the tip than the behavior of 
the shape docs. The picture is the following: As one moves down from the tip, one first finds the linear region where 
branches are born and eventually start to compete with cachothcr. Going further, the nonlinear region appears after 
a transition region. The competition process is not only still taking place there but it is probably in the highest 
point of activity. Not so far, its effects will be easily seen by the observation of some already stopped sidebranches. 
During all this way we have moved from the linear to the nonlinear region, but we are still in the region we called 
A. Further down, the competition process between branches is finished and the surviving ones have no opposition in 
their neighborhood to keep growing as free dendrites. We are then in region B. 

Finally, we have considered both surface tension and kinetic dominated dendrites. Although the different behaviors 
of the studied parameters arc observed when the undercooling is changed, this cannot be associated to the type of 
dendrite. What we can only assure is that the linear and A regions in surface tension dendrites will always be larger 
than in the kinetic ones, but only because of the larger diffusion length and not because of the main mechanism which 
determines them. 

These results offer some new insight into the understanding of a fully developed dendrite, and in particular are 
of relevant importance to distinguish between low and high undercooling dendrites and two-dimensional and three- 
dimensional dendrites. Finally we should remark that it would be of the most great interest to have more experimental 
results available in the high undercooling regime, in particular characterizing the nonlinear regions of the dendrite. 
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